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(N : ABSTRACT 



The non-spherically symmetric transport equations and exact thermal evolu- 
tion model are used to calculate the transient thermal response to pulsars. The 



^ ! three possible ways of energy release originated from glitches, namely the 'shell' 



'ring' and 'spot' cases are compared. The X-ray light curves resulting from the 



O ' thermal response to the glitches are calculated. Only the 'spot' case and the 
O ' 

O ■ 'ring' case are considered because the 'shell' case does not produce significant 

Qh' modulative X-rays. The magnetic field [B) effect, the relativistic light bending 

2 ' effect and the rotational effect on the photons being emitted in a finite region 

■ are considered. Various sets of parameters result in different evolution patterns 



of light curves. We find that this modulated thermal X-ray radiation resulting 
from glitches may provide some useful constraints on glitch models. 



Subject headings: dense matter-stars: evolution-star: interiors-stars: neutron- 
stars: X-rays 



1. Introduction 

Pulsar glitches are believed to be mainly either starquake-driven or superfluid-driven. 
The starquake mechanism was introduced by Ruderman (1969). He suggested that glitches 
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could be caused by gravity. As a pulsar spins down due to dipole radiation, centrifugal force 
on the crust decreases and gravity pulls the crust towards a less oblate shape. Since pulsars 
are believed to be neutron stars with a solid crust, subsequent change in stellar shape induces 
stress in the crust until the maximum yield strength is reached. Sudden relaxation of this 
stress brings the stellar shape to equihbrium. The glitch energy is to be released in a small 
volume at the weak regions in the solid crust. This leads to the 'spot' proposed 
by Van Riper et al. (1991). This localized heating process causes an uneven heating of the 
surface. As the pulsar rotates, the area of thermal photon emission facing the observer is 
changing. Together with the gravitational bending effect, the emission of thermal X-ray due 
to glitches will be at a particular phase in the light curve. Hence a periodic modulation of 
the X-ray light curve will be observed. 

The superfluid-driven glitch mechanism was proposed by Anderson (1975). When the 
pulsar spins down, the vortices become unpinned from the lattice. They corotate with the 
local superfiuid and scatter with the nearby lattices. There are two effects. One is superfiuid 
angular momentum being transferred to the crust so that the crust spins up. The other 
is frictional heating being produced and causing local energy dissipation. Such a process 
occurs in a ring structure at the rotational equator (Alpar, Anderson, Pines and Shaham 
1984). The detailed study of glitch energy depositing in a ring structure has been studied 
by other authors (eg. Bildsten and Epstein 1989; Epstein and Baym 1992; Link and Epstein 
1996; Jones 1998). 

In reality, the 'spot-like' case and the 'ring-like' case may not necessarily be mutually 
exclusive. After crust-breaking, which will cause a 'spot' case, the whole crust will oscillate 
briefly about its new equilibrium conflguration (Baym and Pines 1971). The vortices which 
are originally pinned to the crust will 'feel' this oscillation and become unpinned from the 
crustal lattice. These depinned vortices start to scatter outward and transfer their differential 
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rotational energy to the crust. Therefore, the 'ring-hke' ghtch can follow. As a result, we 
believe that the energy release in a glitch is mainly composed of two components, namely 
the release of elastic energy in a starquake and the differential rotational energy released 
due to the depinning of vortices. On the other hand, several authors (Van Riper et al. 1991; 
Chong and Cheng 1994; Hirano et al. 1997) assume the glitch energy to be released in a 
spherical shell at a certain density inside the pulsar, namely the 'shell' case, though this 
does not seem to be realistic. There are other possible origins of the 'spot' case of energy 
release (e.g. Srinivasan et al. 1990; Ruderman 1991a; Link, Franco and Epstein 1998) and 
the 'ring' case of ghtch energy release (e.g. Bildsten and Epstein 1989; Epstein and Baym 
1992; Jones 1998). 

It has been suggested that the transient X-ray emission resulting from glitches will pro- 
vide a good method to determine the equation of state for pulsars (neutron stars). However, 
the most active glitching pulsars (e.g. Vela, PSR 1706-16) are young neutron stars whose 
interior temperature is high (~ 10® K)., this makes the overall luminosity variation difficult 
to detect (Van Riper et al. 1991; Chong and Cheng 1994). Cheng, Li and Suen (1998) argue 
that if a good fraction of glitch energy is released in a small volume, then instead of heating 
up the entire stellar surface, even a small fraction of glitch energy can heat up a small area 
of the stellar surface drastically. Therefore, although after a glitch the total X-ray intensity 
varies very little, a very distinctive hot spot may be detected. They suggest that by observ- 
ing the transient X-ray pulses, the equations of state of neutron stars can be determined. 
However, they have not considered such important effects as the relativistic light bending 
(Pechenick et al. 1983) and the magnetic field effect (Page 1995) which can significantly 
affect the intensity and the pulse shape of the transient X-rays resulting from glitches. 

In Section ^ we summarize the input physics and the relativistic non-spherical sym- 
metric thermal transport and energy balance equations that are used in determining the 
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cooling process following a glitch. In Section we apply the scheme mentioned in Section ^ 
to the three cases, namely the 'spot', 'ring' and 'shell' cases, and compare the so-obtained 
temperature and luminosity profiles. In Section ^, we discuss the magnetic field effect, the 
relativistic light bending effect and the rotational effect. In Section we calculate the 
expected periodic modification to thermal X-ray pulses that are emitted due to 'spot-like' 
and 'ring-like' glitches resulting from the effects mentioned in Section ^. In Section ^ we 
briefly summarize our results and discuss the detectability of the thermal X-ray profile being 
modified due to pulsar glitches. 



2. Physics inputs 

No matter what the origins of the glitches are, glitches can be simulated by energy 
deposition in particular regions in the pulsar. The excess energy propagates in all directions. 
In this paper, we would like to know what fraction of glitch energy, which gives the transient 
emission of electromagnetic radiation from the stellar surface, can be observed. In calculating 
the thermal afterglow effect of glitches, the properties of a neutron star considered to be the 
major factors affecting the energy flow are the equations of state, the composition and the 
initial temperature profile. 



2.1. Neutron star structure 

The structure of a neutron star can be constructed by the general relativistic hydrostatic 
equilibrium equation (Tolman-Oppenheimer-Volkoff equation) 

dP G{m + 47rr^P/c^) (p + P/c^) 



dr r{r — 2Gm/c^) 



(1) 
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where P and p are the pressure and mass density at radius r respectively, and 

rr 

m — 4:nr^pdr (2) 

^0 

is the gravitational mass inside radius r. G and c are the gravitational constant and the 
speed of light in a vacuum respectively. 

We will only calculate the temperature profiles between the radii with densities between 
10^ gjcvrc' and the nuclear density p^v- Prom 10^ gjcrr? to the stellar surface, we assume that 
the region between these two densities reaches equilibrium quickly, so that the temperatures 
at these two densities are related by using the formula proposed by Gudmundsson et al. 
(1983). In this paper, since different equations of state give rise to different kinds of structure 
for a neutron star, the major equation of state (EOS) used is UT, which is based on the 
combined Hamiltonians consisting of UV14 and TNI models (Largaris and Pandharipande 
1981; Priedman and Pandharipande 1981; Wiringa and Piks 1988). It is a moderately stiff 
equation of state, when compared to the softer BPS model of Baym et al. (1971), and the 
stiffer PPS model of Pandharipande et al. (1976). 



2.2. Composition 

The ion, neutron and proton mass fraction, electron fraction, the mass number and 
proton number in the crust region denoted by Xion, Xp, Yg, A and Z respectively 
are obtained from Lattimer et al. (1985). They are important for calculating the neutrino 
emissivity, heat capacity and thermal conductivity as presented in the following subsections. 
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2.3. Neutrino emissivity 

The most important neutrino emission process for pulsars with age ^ 10^ years is 
electron bremsstrahlung in the crust (Flowers and Itoh 1976, 1979) as 

72 

Qions ^ 2.1 X 10^° — ( — )T96 erg cm'^ (3) 
A pn 

where Tg is the temperature in units of 10^ K and pn — 2.8 x 10^^ g/cm^ is the nuclear 
density. Other neutrino emissivities (cf. Chong and Cheng 1994 for a brief review) are also 
included. 



2.4. Heat capacity (C^) 

In the crust, there are extremely relativistic degenerate electrons, some non-relativistic 
neutrons and ions, but there are no free protons. The scheme for calculating heat capacities 
follows Chong and Cheng (1994) and Cheng, Li and Suen (1998). 



2.5. Thermal conductivity (K) 

Thermal conductivity is fitted according to the work of Itoh et al. (1984) in solid phase 
and that of Itoh et al. (1983) in liquid phase. In regions of p > 1.311 x 10^^ g/cm^, hnear 
extrapolation is used. 



2.6. Initial temperature 

Before a glitch occurs, the temperature of the pulsar should be in equilibrium. For 
young and middle-aged pulsars, in which glitches can occur, the core temperatures are about 
~ 10^ — 10^ K. Since different EOS give rise to different kinds of structure, the equihbrium 



- 7- 



temperature profile of a pulsar depends on the EOS. 



2.7. Non-spherically symmetric general relativistic transport and energy 

balance equations 

According to Cheng et al. (1998), the general relativistic thermal transport equation 
and the energy balance equation are given by 

^(r^F.e^*) + -^^(sin^F,) + -^^F, = -(a^ + e*g.) (4) 



dr r sin 9 89 r sin 9 d(f) dt 



and 



|:M =-'^Fr (5) 
l&(Te') --"^Fe (6) 

where T is the temperature, e* and are the redshift factor and length correction factor 
respectively, and F^, Fq and F^ are the heat fiuxes along r, 9 and (f) directions respectively 



2.8. Heat input 

As we have mentioned in the introduction, there are likely two types of energy release 
during the glitches, namely, the elastic energy of the crust (Baym and Pines 1971) and the 
differential rotation energy between the crustal superfiuid and the solid crust. In the former 
case, the energy is released in a locahzed volume, which is referred to as the 'spot' case. In 
the latter case, the energy is released in the equatorial plane and hence is called the 'ring' 
case. The amount of energy released in these two cases is estimated as follows: 



-8- 



2.8.1. 'Spot' case 

The energy released in a starquake comes from the rehef of strain energy and is estimated 
to be (Baym and Pines 1971; Ruderman 1991b; Cheng et al. 1992) 

where ji [~ 10^^(p/10^^ g/cm^)~'^/^ dynes cm~^ where p is the mass density] is the mean 
shear modulus, Vcmst is the volume of the crust from which energy is released, and 9max 
[~ 10^^ — 10~^ for a pure coulombic lattice and 10"'^ — 10"'^ for an impurity dominated lattice 
(Smolukowshi 1970)] is the maximum strain angle that the crust can withstand without 
cracking. If the glitch occurs at p ~ 10^^ g/cm^, the strain energy released is estimated as 
^Estrain ~ 10^°(6'maa:/10~^)^ cvg. Hence, the estimated AEstrain ^ 10''° erg. 

It has been argued that the magnetic pressure of magnetars, which are neutron stars 
with extremely strong magnetic fields (~ 10^^ G), is strong enough to break the crust during 
the evolution of magnetic field and may encourage glitching with an amplitude ^ ~ 10~^ 
(Thompson and Duncan 1996). In fact, Heyl and Hernquist (1999) have found evidence 
for glitches occuring in possible magnetar candidates, IE1048. 1-5937 and IE2259-I-586, with 
amplitude ^ ~ 10~^. If a large fraction of magnetic energy is dissipated inside the crust, 
the heat dissipation resulting from the glitch can be estimated as ~ 10^° erg usmg 

the typical parameters of magnetars (Q ~ 1 rad/s, ^ ~ 10~^ and /^^ ~ 10^^ g crri^ where 
li, is the total moment of inertia of the star) . 

2.8.2. 'Ring' case 

In a superfiuid-driven glitch originated from a sudden transfer of angular momentum 
from the inner crust superfluids to the crust, the angular momentum loss for the crustal 
superfluid is Icr^^s where Icr is the moment of inertia for the crustal superfluid and Silg is 
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the angular velocity change of crustal superfluid before and after the ghtch. The angular 
momentum loss for the charged component, including the stellar core, which is strongly 
coupled to the solid crust via electron-magnetized vortex scattering (Alpar, Langers and Sauls 
1984), and the solid crust, is /c/iAQ, where AQ is the observed angular momentum jump of 
the glitch and Ich is the moment of inertia for the charged component. By angular momentum 
conservation A J = Ich^^ = hr^^s and the energy dissipated due to the loss of differential 
rotation energy between the crustal superfluid and the charged component is 

AE = ^[/eftO' + ler^l] - ^[Ich{^ + A^f + hri^s - S^sf] = AJ(1], - = AJ^iag. (9) 

where fliag is the angular speed difference between the crust and the superfluid. In the 
nuclear pinning region, Qiag is ~ 1 — 100 rad/s (Alpar, Cheng and Pines 1989). In the 
interstitial pinning region, it is ^ 0.1 rad/s (Link and Epstein 1991). According to the 
postglitch relaxation (Epstein, Van Riper and Link 1992; Alpar, Chau, Cheng and Pines 
1993), Ich ^ Icr, AJ = Ich^^ — {h — /cr)AO ~ I^Afl. For a typical neutron star with a 
giant glitch AQ ~ 10~^ rads~^, the moment of inertia of the star is ~ 10^^ g cm? and the 
estimated energy release in nuclear pinning regions is ~ 10^-^ — 10^^ erg, while that in the 
interstitial pinning regions is ~ lO"^'' erg. 

In general, the starquakc glitches (the 'spot-like' cases) can occur anywhere within the 
crust where a coulomb lattice exists, whereas superfluid- driven glitches (the 'ring' cases) 
can only occur in the inner crust where neutron superfluid and coulomb lattices coexist. 
Therefore, they must occur in the inner crust at densities between 10^^ and 2 x 10^^ gjcrr?. 
For simplicity, the 'ring' case is assumed to occur at Pgutch ~ 

1013 _ 1014 g/cm^. Other 

cases can occur at Pgutch ~ 10^^ — 10^^ gjcrr? . The amount of energy liberated, AE, is thus 
estimated to be between ~ 10^*^ erg and 10^^ erg. 
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3. Comparison of 'shell', 'ring' and 'spot' cases 

A finite difference method is used to solve the equations of the previous section. For 
the 'shell' case, spherical symmetry can be confidently assumed. As for the 'spot' case, the 
cell at which energy release takes place can be treated as a 'pole'. The direction joining the 
centre of the pulsar and the 'pole' is defined to be the z-direction. Since the rotation of the 
pulsar is slow, the direction of rotational axis can be neglected and azimuthal symmetry can 
be assumed. For the 'ring' case, the angular velocity direction is taken to be the i'-direction, 
so azimuthal symmetry can be assumed again. Therefore, a 2-dimensional grid with A^^ x Nq 
cells is used, where A^^ and Nq are, for computational convenience, taken to be 100 and 90 
respectively. 

For the 'spot' case, the release of energy is in one cell of the grid. For the 'shell' case, 
the release of energy is in a layer of the grid. Whereas for the 'ring' case, the release of 
energy is around the equator with a height of 4°, i.e. ~ 0.69 km for a UT-star, which is 
roughly the thickness of the crust. 

The main parameter used for comparison of these various cases, is the total surface 
luminosity as observed by an observer at infinity, i.e. Lf. As L oc T^, the increase in 
temperature is reflected in the luminosity profile. 

The temperature profile for a 'ring' case glitch is shown in Fig. The distribution of 
surface temperature with respect to the angle from the rotational axis is shown in Fig. ^ 
The distribution of surface temperature with respect to the angle from the spot for the 'spot' 
case is shown in Fig. ^. The evolution for the 'ring' case is faster than the 'spot' 
presented in Cheng, Li and Suen (1998), but it is slower than the 'shell' case as presented in 
Chong and Cheng (1994). 

Since the heat propagation time is proportional to the heat capacity, which is approx- 
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imately proportional to the temperature, i.e. oc T, and inversely proportional to the 
thermal conductivity, which is approximately inversely proportional to the temperature, i.e. 
K (X ^, the time that the glitch energy requires to reach the stellar surface is roughly propor- 
tional to the temperature squared, i.e. t oc ^ oc (Hirano et al. 1997). It can be seen that 
the 'shell' case gives the fastest response to the glitch because the same amount of energy 
is deposited in a larger volume, the localized heating effect is averaged out. Therefore, the 
temperature in each particular area is not so high and the heat energy can reach the surface 
faster than in the other two cases. 

The duration of the pulse on the surface for the 'ring' case is also between that for the 
'shell' case and the 'spot' case. This is similar to Chong and Cheng (1994) 's 'shell' case, with 
a larger amount of energy deposited. The pulse will last a longer time because more heat 
is deposited in a particular volume and this will increase the diffusion time required to the 
next cell. Therefore, the pulse due to a 'spot-like' glitch will last longer, while the second 
shortest one is the 'ring' case and lastly the 'shell' case. 

On the other hand, the peak luminosity for the 'spot' case is the highest because the 
localized heating effect is transferred to the surface. Also, the increase in luminosity lasts the 
longest time for the 'spot' case because the heat diffusion speed is the lowest (c.f. Fig. Ijand 
Fig. 1^). The increase in luminosity is smaller for a cooler core than for a hotter core. The 
heat content of a pulsar is significantly enhanced by a glitch if the energy liberated is large 
compared to the original heat content (Chong and Cheng 1994). This can be understood by 
the fact that luminosity is proportional to T^. Hence the temperature difference is magnified. 
However, the fractional increase ^^'^'^°^)-^!'^^) of a cooler core is ~ 1.3 and only ~ 0.014 for 
a hotter core. Therefore, it is much easier to observe the transient X-ray pulse from a cooler 
star than from a hotter star. 

In Fig. H, we can see that a softer EOS like BPS gives rise to a faster response and a 
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higher peak. This is because a softer EOS gives a thinner crust. Heat diffuses to the surface 
faster because the distance to the surface is smaller. The deposition of energy to the core 
will also be smaller and a larger amount of the released energy is transported to the surface. 

4. Thermal X-ray modification 

Andersen and Ogelman (1997) have already proposed that thermal afterglow caused by 
transient energy releases in a neutron star can alter its X-ray pulse shape by heating a portion 
of the crust so that more thermal X-rays are emitted at a particular phase. By comparing 
such kinds of changes in X-ray pulses with the model results, constraints on glitch models 
can hopefully be made. In the previous section, the temperature profiles on the surface of 
the pulsar for various models are obtained. However, in order to calculate thermal X-ray 
light curves resulting from glitches the gravitational lensing effect and the magnetic field 
effect must be considered. 

4.1. The effect of magnetic field 

It is well-known that the thermal conductivity in a magnetized neutron star depends 
on the angle between the heat current and the magnetic field because both the heat con- 
duction and the opacity coefficient depend on the magnetic field direction (Hernquist 1984; 
Yakovlev and Urpin 1980; Yokovlev 1982; Tsuruta 1986; Schaaf 1987,1988, 1990). The 
anisotropic conductivity in the crust creates a distribution of the temperature over the stellar 
surface. Fortunately, it has been shown that an approximate solution to the thermal diffu- 
sion equation that accounts for this magnetic field effect can be reduced to a one-dimensional 
problem and the surface temperature distribution is given by (e.g. Page 1995). 

rj;„„,(S, T,, e^) = T^,JB = 0, Ob = 0)(cos2 Ob + xl{B, T,) sin^ 6^) (10) 
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with Xo = K±/K^i, the ratio of the thermal conductivities perpendicular and parallel to the 
magnetic field, which according to Greenstein and Hartke (1983), is assumed to be constant 
within the envelope, xo is further established as "■""^ ^ 

Page (1995). Qb is the 

angle between the local field and the radial direction. Here, only the dipolar magnetic field 
case is considered. Ti, is the temperature at a density of 10^° g/cw?, Tg ani is the surface 
temperature of the anisotropic case and T^^jso is the surface temperature of the isotropic (no 
magnetic field) case. 

4.2. The effect of relativistic light bending 

Harding and Muslimov (1998) modeled the soft thermal X-ray profile by noticing that 
the magnetic polar caps of a cooling neutron star are slightly hotter than the rest of the 
stellar surface because of the strong magnetic field effect on heat transport in the surface 
layers. A similar approach is used to analyze the 'spot' case in this paper, but the origin of 
the excessive heating is different. In fact, in their case, a steady X-ray pulse is expected but 
in our case, a transient X-ray pulse occurs. In this paper, the effect of the magnetic field is 
ignored. 

The following scheme of calculating the X-ray light curves resulting from the gravita- 
tional bending and the stellar rotation is adapted from Pechenick et al. (1983). Various 
model temperature profiles of the 'spot' case are incorporated to determine the energy fiux 
emitted from different positions of the pulsar. Since a neutron star has a large mass and a 
small radius, it is necessary to consider the gravitational deflection of the emitted photons 
(c.f. Fig. When the photons are emitted at an angle 5 from the observer's direction, it will 
seem to the observer that they are emitted at an angle 9' . The relationship between 5 and 
9' for different EOS are plotted in Fig. ||. The 'flat' curve, which ignores the gravitational 
bending effect, is used for comparison. A softer EOS gives rise to a stronger lensing effect. 
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This relation is generated from the equations given by Pechenick et al. (1983): 

i + p = e' (11) 

GM 2 

/5 = //^'[0"-(l-2^)-']"''^^ (12) 



. _i Dsin^ 
= sm 



R \ 1- 



2GM 



2GM 



) (13) 



Dc2 

where D is the distance to the observer and the angles are defined as in Fig. |^. 

To avoid the complex procedures in solving the above three equations, some approx- 
imations can be made. Since is extremely small, and ^ is also very small, 5 can 

GM 2 2 

be treated as ^(1 - ^) and {3 ^g^{^ + T^ gm/c^ )^) which is approximately 

f + ^(g^)^ = f + I(f As a result, 9' ^ P ^ 5^+ |(f )^. Hence, 9' -5^ f|5. 
Actually, this is a good approximation for 6 < 60°. 

Since the rotation speed of a pulsar is much less than the speed of light, it can be treated 
clS db slowly rotating rigid body, in which case 

cos(^^o) = sin(7) sin(7o) cos{Qt) + cos(7) cos(7o) (14) 

where 9q is the angle between the centre of 'spot' region and the observer, 7 and 70 are the 

position angles from the z-axis to the centre of 'spot' region and the observer respectively 

as described in Fig. ^. Qt is the azimuthal angle at time t and Q is the stellar angular 

velocity. So far, we have considered the emission region as a point, however, in reality, the 

emission region is of a finite area. If the photon emission occurs at a surface area defined in 

a cone with an angle a at the stellar surface, the result is then generated from the following 

equations according to Pechenick et al. (1983): 

_i cosa-cos^ocos^ 

h[9; a, 9o] =2 cos ■ n ■ n 15 

sm 9o sm 9 

for 9 between 9q ± a. 

9{x) = / - (1 - 2?i)n2]-V2 (16) 
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^0 is defined in the same way as in Eq. |T^. The relative brightness is 



2GM o GM o r^'™- 
A{eo;f,M/R,a) = f{S{x))h[x;a,eo]xdx (17) 

where f{5) = 1 for isotropic emission, f{6) = cos 6 for enhanced emission and f(S) = sin 5 
for suppressed emission, 6 = sin"^(a;/a;maa;)- As ^ is small, sin^ is approximated to ^, so 
X = D^c^/GM and 

Xma. = {m/GM){l - 2GM/R^)-^'^ (18) 

where the angles are as defined in Fig. ^ 

The relationship between A and VLt is plotted in Fig. |T^ for an isotropic emission. A 
softer EOS gives a stronger lensing effect so that when the 'spot' is nearly at the back 
(~ 170°), it can still be observed even though the relative brightness is low. For a UT- 
star, the bending effect can extend the observed angle to nearly 30°. With a larger cone of 
emission, the relative brightness is larger when the 'spot' is facing the observer, but it does 
not have any effect on the maximum angle of deflection. 

According to Pechenick et a/. (1983), the total energy flux observed is: 

Fx = Y: /o(§)'^(^o; /, M/R, a) (19) 

where /q is the energy flux at the surface for different angles from the spin axis and includes 
the factors cos 5 and X); which denotes the summation of the contributions from each cell. 

Including the gravitational bending, together with the slowly-rotating and the flnite 
area effects, the evolution of total brightness for a UT-star with = 10* K, AE = 10^° erg 
and pgiitch = 10^^ g/cm? is plotted in Fig. [TT|. Other sets of parameters give similar shape 
light curves, except that the number of days to reach the maximum will change, as will the 
fractional increase in energy flux. 



- 16 - 



5. Numerical results 

Consider the fractions fp = ^""p — and fx = — which measure the 

nm 'inn ^ 



min 



visibility of the X-ray pulses due to a 'spot-like' glitch in one particular rotation around the 
peak of the luminosity-time graphs. The numerical values are given in Table |l], where the 
magnetic field effect has not yet been included. 

If fp > 0.05, it should be large enough for observation, then the cases with = 10^ K, 
AE = 1039 erg and Pgutch = W g/cm^] = 10^ K, AE = 10^9 erg and Pgutch = 
10^° g/cm^; = 10^ K, AE = 10^° erg and Pgiuch = 10" g/cm^; and = 10^ K, 
AE = 10^^ erg and Pgutch = 10^'^ g/cm? are possible to be observed. 



In the legend of Fig. |T2|, the first number is the angle of 'spot' from the rotational axis 
and the second number is the angle of observer from the rotational axis. Other numerical 
results can be found in Tang (1999). 

Taking the magnetic field effect into consideration, the background light curve is no 
longer a straight line. That is to say, the two magnetic poles have already produced modu- 
lation to the flux profiles. For easy calculation, the background light curves in the figures in 
this second part of the section are generated by using a 10° x 10° grid; the rotational axis Vt 
is taken to be the i'-axis; the magnetic moment p, is taken to be at an angle of 45° from fi; 
the 'spot' case glitch has its energy released at an angle of 60° from VL and 90° azimuthally 
from p. In other words, if p) is the z — y plane, then the spot is in the z — x plane. The 
'ring' case glitch has its energy released along the rotational equator. The observer is at an 
angle of 45° from VL. 

For a 'spot-like' glitch, the glitch effect alone can produce modulation to the X-ray 
profiles as mentioned in the beginning of this section. It can significantly increase the fraction 
fp with some particular sets of parameters. If there is a magnetic field effect, it can also 
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cause a phase shift of the whole hght curve as seen in Fig. 13 to Fig. T7. The results are 
summarized in Table 0. 

For the 'ring' case, if there is no magnetic field, there will be no modulation to the light 
curve due to azimuthal symmetry. Only the observed total flux is increased. If a magnetic 
field effect is also considered, the glitch effect causes the peak of the X-ray profile to rise and 
the pulse shape to shift. However, it depresses the contrast fp- This is easily observed from 
the values in Table ^ and in Fig. |18| to Fig. |19[ 

The rise of the peak flux value can be more drastic for the 'spot' case than the 'ring' 
case with the same set of parameters, especially at low temperature, e.g. Tcore = 10^ K, and 
with a large amount of energy deposition, e.g. AE = 10^"^ erg. The amount of phase shift 
of the light curve is determined by the position at which the glitch energy is released (its 
azimuthal angular distance from fl). For the 'spot' case, since we choose the 'spot' to be 90° 
from fl azimuthally, there is a 90° phase shift in the X-ray light curve when the glitch energy 
reaches the surface (c.f. Fig. |T^). For the 'ring' case, the additional peak results from the 
glitch (c.f. Fig. 1^). For both the 'spot' and the 'ring' cases, the values of fp are greater 
than 0.05. As a result, we expect that the change of light curve due to the glitch should be 
observable. In fact, the change in pulse shapes point to the corresponding sets of parameters. 

It is important to note that with the magnetic field effect, the heat transport becomes a 
3-dimensional problem even for the 'ring' case. The energy released in a glitch can have two 
routes. Part of it is transferred to the core (c.f. Fig. [^). The core is able to quickly resume 
its isothermality and the energy is reradiated isotropically to both directions of the magnetic 
pole and the magnetic equator. According to Harding and Muslimov (1998), the transverse 
heat conductivity is suppressed due to magnetization of electrons, which results in a surface 
temperature at the magnetic pole higher than at the equator. The other part of the glitch 
energy is transferred to the crust, since the heat propagation time is proportional to the 
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temperature squared as mentioned in Section |^. The distance from and the temperature 
of the heated area are the main elements affecting the angular distribution of the thermal 
energy on the surface. 

When the 'ring' case is compared to the 'spot' case, the same amount of energy is 
deposited in a larger volume and the localized heating effect is averaged out. Therefore, 
the temperature in each particular volume element is not so high. Hence, the heat energy 
can reach the surface faster. Comparing Fig. which describes the thermal evolution of a 
'ring-like' glitch, to Fig. 1 in Cheng and Li (1997), which describes the thermal evolution of 
a 'spot-like' glitch, this effect can be easily seen. That means, for a 'ring-like' glitch, a larger 
part of the glitch energy is transferred to the crust than to the core. Since the magnetic 
pole and the magnetic equator are at the same distance from the 'ring' (with the magnetic 
moment fl is at an angle of 45° from Q and the 'ring' is at 90° from Q), and the magnetic 
pole is hotter, energy transport in the magnetic equator direction is faster hence decreasing 
the contrast fp as shown in Table ^. For the 'spot' case, however, a larger fraction of the 
glitch energy is transported to the core, which explains the increase in the contrast fp for 
the 'spot' case. 

For the case of low temperature, Tcore = 10^ K, and large energy deposition, AE = 
Iq42 ^^g^ since the localized heating is enormous, the effect of polar cap heating can be 
viewed simply as a modification to the uniform temperature profile. The magnetic effect is 
absolutely overwhelmed by the glitch effect. However, for a young star like Vela, the core 
temperature is still high (~ 10® K), the rate of energy transport is slower, and hence a 
larger fraction of the glitch energy is transferred to the core. Therefore, the glitch effect is 
not significant both in terms of phase shift and the rise of the peak of flux. 

We define a quantity = /q^*'"'"' (L(t) — LBG)dt, where tmax is the time at which L{t) 
reaches its peak value after the glitch and Lbg is the luminosity before the glitch. This 



- 19 - 



parameter indicates how much ghtch energy is released as a pulse. From the " E^" column 
in Table ^ and Table ^ it can be easily seen that the faster the luminosity (L) reaches its 
peak value, the smaller is the fraction of the glitch energy eventually emitted as photons at 
the surface within 2tmax- We wish to point out that eventually all the glitch energy will be 
radiated through the surface in a time scale much longer than tmax- 



6. Conclusion and discussion 

In this paper, we have employed the relativistic heat diffusion equations to calculate 
the thermal response of a pulsar to a glitch with energy deposited in a ring-like area around 
the rotational equator in the inner crust, in a spherical shell and in a small volume. The 
'ring' case is always the middle one in terms of the response time, the response period and 
the intensity of the response. This is not unexpected as the volume of energy deposition is 
smaller than the 'shell' case but larger than the 'spot' case. Moreover, we have calculated the 
expected light curves for the thermal X-ray regime as observed after a glitch for the 'spot' 
case and the 'ring' case. The 'shell' case is ignored because we believe that the modulation 
to the light curves in this case will be too insignificant to be observed. After a glitch, the 
energy flux emitted from the surface is not uniform. Excess thermal X-rays are emitted at a 
particular phase. This alters the X-ray pulse shape of the pulsar. These curves can be used 
to compare with the actual data obtained to put constraints on the true glitch model as well 
as the interior structure of a neutron star. 

For the prompt afterglow of a glitch to be detectable, the energy released should be 
enormous so that it should re-radiate over a reasonable time scale and the change in intensity 
of the radiation can be observed. Thermal emissions for AE ^ 10^^ erg lie in the soft X- 
ray region (Chong and Cheng 1994). The thermal transients produced by giant glitches in 
nearby pulsars may be observable by Chandra, XMM and Astro-E. These three missions 
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have complementary capabilities in measuring soft X-ray pulses from neutron stars. 

For observations taken following a glitch that triggers thermal afterglow, one should 
be able to detect significant changes in the X-ray pulse shape if the energy release is 
anisotropic. The results will be valuable in refining neutron star equations of state and 
gaining a better understanding of the physics behind neutron star transient energy releases 
(Anderson and Ogelman 1997). 

The defect of this method of finding the internal properties of a neutron star is that if 
the glitch events occur too frequently, it may result in a pile-up of the pulses as well as a 
long-term variation of the total thermal radiation which eventually reduces the detectability 
of the thermal afterglow of the glitch (Li 1997). In the calculations in this paper, several 
assumptions have been made: the interstellar absorption, the magnetospheric effects and 
the possibility of other diffusion mechanisms besides conduction in the pulsar have been 
neglected; the structure of the pulsar is assumed to be inert to the thermal change of the 
pulsar; the thermal X-ray emitted from the surface is assumed to be blackbody radiation 
alone. In future study, we shall take these factors into account. 

We thank J.S. Heyl, D. Page and S. Tsuruta for useful comments and discussion and 
Terry Boyce for a critical reading. This work is partially supported by a RGC grant from 
the Hong Kong Government and a Croucher Foundation Senior Research Fellowship. 
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Table 1: The fractions fp and fx are calculated for various sets of parameters with the 'spot' 
case of UT equation of state with no B field and spot and the observer both lying along the 
rotational equator. 



Tcore (-^) 


AE (erg) 


Pgiitch {g/cm^) 


fp 


/t 




1038 




0.0002 


0.0395 


10« 


1039 




0.0033 


0.3400 


10^ 


1040 




0.0517 


1.4493 


10« 


10^1 




1.8673 


4.9635 


10^ 


10^2 




50.8427 


12.6201 


10« 


1038 


10^° 


0.0033 


0.3411 


10« 


1039 


10^° 


0.0814 


1.7365 


10« 


1040 


3 X 10^^ 


0.0097 


0.6568 


108 


10^2 


10^=^ 


0.0776 


1.7047 


10^ 


1038 


10^1 


0.0183 


0.9121 


10^ 


1039 


10^1 


0.3670 


3.8872 


10^ 


10^2 


10^2 


66.2741 


13.5537 


10^ 


10^2 


10^3 


9.0100 


7.8376 
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Table 2: The fractions fp are calculated for various sets of parameters for the 'spot' case of 

— * — * — * 

UT equation of state with B at 45° from Q, the 'spot' located at an angle of 60° from Q. and 
90° azimuthally from B, and the observer at an angle 45° with Q. The 'new peak' represents 
the maximum flux value after the glitch and the 'old peak' represents the value before the 
glitch. The value Ej — J^*"""" {L(t) — LBG)dt is the amount of energy that is emitted as 
photons at the surface, tmax is the time at which the luminosity has its peak value after the 
glitch. Lbg is the background luminosity. 
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0.49 
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10^5 


0.40 


0.49 
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Table 3: The fractions fp are calculated for various sets of parameters for the 'ring' case of 

— * 

UT equation of state with jl as well as the observer making an angle 45" with Q. The 'new 
peak' represents the maximum flux value after the glitch and the 'old peak' represents this 
value before the glitch. The value Ej = /q^*'"""" — LBG)dt is the amount of energy that 
is emitted as photons at the surface, tmax is the time at which the luminosity has its peak 
value after the glitch. Lbg is the background luminosity. 
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Figure Captions 

Fig. 1. — The redshifted temperature profile of a UT-star witli Tc = 10^ K, AE = 10*^^ erg 
at p glitch = 10^^ g/cm^ for the 'ring' case with respect to the radius from the centre of the 
star. The top left panel indicates the number of days after the ghtch. 

Fig. 2.— The redshifted surface temperature of a UT-star with Tc = 10^ K, AE = 10^^ erg 
at Pgiitch — 10^^ g/cm^ for the 'ring' case with respect to the angle from the rotational axis. 

Fig. 3.— The redshifted surface temperature of a UT-star with = 10^ K, AE = 10^^ 

at Pgiitch = 10^'^ g/cnv^ for the 'spot' case with respect to the angle from the rotational axis. 
The numbers besides the lines are the number of days after the glitch. 

Fig. 4. — The luminosity of a UT-star with Tc = 10^ K, AE = 10*^^ erg at Pghtch = 
10^^ g/cm^. The three cases-'shell', 'ring' and 'spot'-are compared. 

Fig. 5.— The luminosity of a UT-star with Tc = 10^ K, AE = 10^^ p^^.^^^ = 

10^^ g/cm^. The three cases-'shell', 'ring' and 'spot'-are compared. 

Fig. 6. — The luminosity profiles of a BPS-star, a UT-star and a PPS-star with Tc — 10^ K, 
AE — 10^^ erg at Pgutch — lO-*^^ g/cm^ for the 'ring' case. 

Fig. 7. — Schematic illustration of gravitational lensing effect. Angles are defined as shown. 
Fig. 8. — Gravitational lensing effect for different EOS. 
Fig. 9. — Angles used in calculating the rotational effect. 

Fig. 10. — Finite area for the 'spot' is considered for different EOS: BPS, PPS and UT. The 
angular radius of the emission cone is a = 5°. 

Fig. 11. — The evolution of the total brightness for a UT-star with Tc — 10^ K at Pgiitch — 
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10^^ gjcrr? with ISE — 10^*^ erg for the 'spot' case. 

Fig. 12.— The total brightness for a UT-star with = 10^ K at pgutch = 10^^ gjcn? with 
I\E = 10^2 erg 418 days after a 'spot-hke' ghtch. 

Fig. 13. — The evolution of the total brightness for a UT-star with — 
10^^ gjcrr? with l\E — 10^° erg for the 'spot' case with magnetic field B 

Fig. 14. — The evolution of the total brightness for a UT-star with Tc = 

10^^ g/crrt" with /\E = 10^^ erg for the 'spot' case with magnetic field B 

Fig. 15. — The evolution of the total brightness for a UT-star with = 
10^^ g/cm^ with AE = 10"^° erg for the 'spot' case with magnetic field B 

Fig. 16. — The evolution of the total brightness for a UT-star with Tc = 
10^^ g/cm^ with /\E — 10^^ erg for the 'spot' case with magnetic field B 

Fig. 17. — The evolution of the total brightness for a UT-star with = 
10^^ g/cm^ with AE — 10^^ erg for the 'spot' case with magnetic field B 

Fig. 18. — The evolution of the total brightness for a UT-star with Tc — 
10^2 g/cm^ with AE = 10"^° erg for the 'ring' case with magnetic field S = 3 x 10^^ q 

Fig. 19. — The evolution of the total brightness for a UT-star with Tc — 10^ K at Pgutch — 
10^2 g/cm^ with AE = lO''^ erg for the 'ring' case with magnetic field S = 3 x 10^^ q 

Fig. 20. — The evolution of the total brightness for a UT-star with Tc — 10^ K at Pgutch — 
10^^ gjcrr? with AE — 10^^ erg for the 'ring' case with magnetic field S = 3 x 10^^ G. 



10^ X at p glitch = 
= 3 X 10^2 Q 

10^ X at p glitch = 
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10^ K at p glitch = 
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10^ K at p glitch = 
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